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The theory of nuclear forces has made great progress since the turn of the millenium using 
the framework of chiral effective field theory (ChEFT). The advantage of this approach, 
which was originally proposed by Weinberg, is that it has a firm basis in quantum- 
chromodynamics and allows for quantitative calculations. Moreover, this theory generates 
two-nucleon forces (2NF) and many-body forces on an equal footing and provides an 
explanation for the empirically known fact that 2NF 3> 3NF ^ 4NF. I will present 
the recent advances in more detail and put them into historical context. In addition, I 
will also provide a critical evaluation of the progress made including a discussion of the 
limitations of the ChEFT approach. 
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1. Introduction and historical perspective 

The theory of nuclear forces has a long history (cf. Table 1). Based upon the seminal 
idea by Yukawa/ first field-theoretic attempts to derive the nucleon-nucleon (NN) 
interaction focused on pion-exchange. While the one-pion exchange turned out to 
be very useful in explaining NN scattering data and the properties of the deuteron, 
multi-pion exchange was beset with serious ambiguities. Thus, the "pion theories" 
of the 1950s are generally judged as failures — for reasons we understand today: pion 
dynamics is constrained by chiral symmetry, a crucial point that was unknown in 
the 1950s. 

Historically, the experimental discovery of heavy mesons in the early 1960s saved 
the situation. The one-boson-exchange (OBE) modeP emerged which is still the 
most economical and quantitative phenomenology for describing the nuclear force. '^■^ 
The weak point of this model, however, is the scalar-isoscalar "sigma" or "epsilon" 
boson, for which the empirical evidence remains controversial. Since this boson is 
associated with the correlated (or resonant) exchange of two pions, a vast theoretical 
effort that occupied more than a decade was launched to derive the 27r-exchange 
contribution to the nuclear force, which creates the intermediate range attraction. 
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Table 1. Seven Decades of Struggle: The Theory of Nuclear Forces 



1935 


Yukawaj IVleson Theory 




The "Pion Theories" 


1950's 


One-Pion Exchange: o.k. 




Multi-Pion Exchange: disaster 




Many pions = multi-pion resonances: 


1960's 


<T, p, LJ, ... 




The One-Boson-Exchange Model 




Refine meson theory: 


1970's 


Sophisticated 27r exchange models 




(Stony Brook, Paris, Bonn) 




Nuclear physicists discover 


1980's 


QCD 




Quark Cluster Models 




Nuclear physicists discover EFT 


1990's 


Weinberg, van Kolck 


and beyond 


Back to Meson Theory! 




But, with Chiral Symmetry 



For this, dispersion theory as well as field theory were invoked producing the Paris^ 
and the Bonn-^'^ potentials. 

The nuclear force problem appeared to be solved; however, with the discovery of 
quantum chromo-dynamics (QCD), all "meson theories" were relegated to models 
and the attempts to derive the nuclear force started all over again. 

The problem with a derivation from QCD is that this theory is non-perturbative 
in the low-energy regime characteristic of nuclear physics, which makes direct solu- 
tions impossible. Therefore, during the first round of new attempts, QCD-inspircd 
quark models'^ became popular. These models are able to reproduce qualitatively 
and, in some cases, semi-quantitatively the gross features of the nuclear force. 
However, on a critical note, it has been pointed out that these quark-based ap- 
proaches are nothing but another set of models and, thus, do not represent any 
fundamental progress. Equally well, one may then stay with the simpler and much 
more quantitative meson models. 

A major breakthrough occurred when the concept of an effective field theory 
(EFT) was introduced and applied to low-energy QCD.^*^ 

Note that the QCD Lagrangian for massless up and down quarks is chirally 
symmetric, i. c., it is invariant under global flavor SU{2)l x S'C/(2)/j equivalent to 
SU{2)v X SU{2)a (vector and axial vector) transformations. The axial symmetry is 
spontaneously broken as evidenced in the absence of parity doublets in the low-mass 
hadron spectrum. This implies the existence of three massless Goldstone bosons 
which arc identified with the three pions (7r^,7r°). The non-zero, but small, pion 
mass is a consequence of the fact that the up and down quark masses are not exactly 
zero either (some small, but explicit symmetry breaking). Thus, we arrive at a low- 
energy scenario that consists of pions and nucleons interacting via a force governed 
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by spontaneously broken approximate chiral symmetry. 

To create an effective field theory describing this scenario, one has to write down 
the most general Lagrangian consistent with the assumed symmetry principles, par- 
ticularly the (broken) chiral symmetry of QCD.^° At low energy, the effective degrees 
of freedom are pions and nucleons rather than quarks and gluons; heavy mesons and 
nucleon resonances are "integrated out". So, the circle of history is closing and we 
are back to Yukawa's meson theory, except that we have learned to add one impor- 
tant refinement to the theory: broken chiral symmetry is a crucial constraint that 
generates and controls the dynamics and establishes a clear connection with the 
underlying theory, QCD. 

It is the purpose of the remainder of this contribution to describe the EFT 
approach to nuclear forces in more detail. 

2. Chiral perturbation theory and the hierarchy of nuclear forces 

The chiral effective Lagrangian is given by an infinite series of terms with in- 
creasing number of derivatives and/or nucleon fields, with the dependence of each 
term on the pion field prescribed by the rules of broken chiral symmetry. Applying 
this Lagrangian to NN scattering generates an unlimited number of Feynman dia- 
grams. However, Weinberg showed^^ that a systematic expansion exists in terms of 
{Q/A^y, where Q denotes a momentum or pion mass, 1 GeV is the chiral 

symmetry breaking scale, and > (cf. Fig. 1). This has become known as chi- 
ral perturbation theory (xPT). For a given order v, the number of terms is finite 
and calculable; these terms are uniquely defined and the prediction at each order is 
model-independent. By going to higher orders, the amplitude can be calculated to 
any desired accuracy. 

Following the first initiative by Weinberg, pioneering work was performed by 
Ordofiez, Ray, and van Kolck^^'^^ who constructed a NN potential in coordinate 
space based upon xPT at next-to- next- to-leading order (NNLO; v = 3). The results 
were encouraging and many researchers became attracted to the new field. Kaiser, 
Brockmann, and Weise^^ presented the first model-independent prediction for the 
NN amplitudes of peripheral partial waves at NNLO. Epelbaum et al}^ developed 
the first momentum-space NN potential at NNLO, and Entem and Machleidt""^^ 
presented the first potential at N^LO (i^ = 4). 

In xPT, the NN amplitude is uniquely determined by two classes of contribu- 
tions: contact terms and pion-exchange diagrams. There are two contacts of order 

[0((5")] represented by the foTir-nuclcon graph with a small-dot vertex shown in 
the first row of Fig. 1. The corresponding graph in the second row, four nucleon legs 
and a solid square, represent the seven contact terms of 0{Q'^). Finally, at 0{Q'^), 
we have 15 contact contributions represented by a four- nucleon graph with a solid 
diamond. 

Now, turning to the pion contributions: At leading order [LO, 0{Q^), v = 0], 
there is only the wellknown static one-pion exchange, second diagram in the first row 
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2N Force 3N Force 4N Force 




Fig. 1. Hierarchy of nuclear forces in xPT. Solid lines represent nucleons and dashed lines pions. 
Further explanations are given in the text. 

of Fig. 1. Two-pion exchange (TPE) starts at next-to-leading order (NLO, i/ — 2) 
and all diagrams of this leading-order two-pion exchange are shown. Further TPE 
contributions occur in any higher order. Of this sub-leading TPE, we show only two 
representative diagrams at NNLO and three diagrams at N^LO. The TPE at N^LO 
has been calculated first by Kaiser. All 2tt exchange diagrams/contributions up 
to N'^LO are summarized in a pedagogical and systematic fashion in Ref. 18 where 
the model-independent results for NN scattering in peripheral partial waves are also 
shown. 

Finally, there is also three-pion exchange, which shows up for the first time at 
N^LO (two loops; one representative Stt diagram is included in Fig. 1). In Rcf. 19 
it was demonstrated that the Stt contribution at this order is negligible. 

One important advantage of xPT is that it makes specific predictions also for 
many-body forces. For a given order of xPT, two-nucleon forces (2NF), three- 
nucleon forces (3NF), ...are generated on the same footing (cf. Fig. 1). At LO, 
there are no 3NF, and at next-to-leading order (NLO), all 3NF terms cancel. 
However, at NNLO and higher orders, well-defined, nonvanishing 3NF occur. 
Since 3NF show up for the first time at NNLO, they are weak. Four-nucleon forces 
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Table 2. /datum for the reproduction of the 1999 np database by famihes of np potentials at 
NLO and NNLO constructed by the Bochum/Juelich group. 





# of np 


Bochum/ Juelich 




Bin (MeV) 


data 


NLO 


NNLO 


0-100 


1058 


4-5 


1.4-1.9 


100-190 


501 


77-121 


12-32 


190-290 


843 


140-220 


25-69 


0-290 


2402 


67-105 


12-27 



(4NF) occur first at N^LO and, therefore, they are even weaker. 
3. Chiral NN potentials 

The two-nucleon system is non-perturbative as evidenced by the presence of shal- 
low bound states and large scattering lengths. Weinberg^^ showed that the strong 
enhancement of the scattering amplitude arises from purely nucleonic intermedi- 
ate states. He therefore suggested to use perturbation theory to calculate the NN 
potential (Fig. 1) and to apply this potential in a scattering equation (Lippmann- 
Schwingcr or Schrodinger equation) to obtain the NN amplitude. We follow this 
philosophy. 

Chiral perturbation theory is a low-momentum expansion. It is valid only for 
momenta Q ^ A-^ ^ 1 GeV. Therefore, when a potential is constructed, all ex- 
pressions (contacts and irreducible pion exchanges) are multiplied with a regulator 
function, 



exp 



(- 
VA 



(1) 



where p and p' denote, respectively, the magnitudes of the initial and final nucleon 
momenta in the center-of-mass system (CMS); and A <^ A-^. The exponent 2n is to 
be chosen such that the regulator generates powers which are beyond the order at 
which the calculation is conducted. 

To what order in xPT do we have to go for sufficient accuracy? To discuss this 
issue on firm grounds, I show in Table 2 the /datum for the fit of the world np 
data below 290 McV for a family of np potentials at NLO and NNLO. The NLO 
potentials produce the horrendous x^/datum between 67 and 105, and the NNLO 
are between 12 and 27. The rate of improvement from one order to the other is very 
impressive, but the quality of the reproduction of the np data at NLO and NNLO 
is obviously totally insufficient for reliable predictions. 

Based upon these facts, it has been pointed out in 2002 by Entem and Mach- 
leidt^^'^^ that NNLO is insufficient and one has to proceed to N^LO. Consequently, 
the first N'^LO potential was created in 2003, which showed that at this order a 
X^ /datum comparable to the high-precision Argonne Vis potential can, indeed, be 
achieved, see Tables 3 and 4. This "Idaho" N^LO potential^^ produces a x^/datum 
= 1.1 for the world np data below 290 MeV which compares well with the x^ /datum 
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Table 3. x'^ /datum for the reproduction of the 1999 np database by various np potentials. 
Numbers in parentheses denote cutoff parameters in units of MeV. 







Idaho 


Bochum/ Juelich 


Argonne 


Bin (MeV) 


# of np 










data 


(500-600) 


(600/700-450/500) 




0-100 


1058 


1.0-1.1 


1.0-1.1 


0.95 


100-190 


501 


1.1-1.2 


1.3-1.8 


1.10 


190-290 


843 


1.2-1.4 


2.8-20.0 


1.11 


0-290 


2402 


1.1-1.3 


1.7-7.9 


1.04 



= 1.04 by the Argonne potential (Table 3). In 2005, also the Bochum/ Juelich group 
produced several N'^LO NN potentials, the best of which fits the np data with 
a x^/datum =1.7 and the worse with a x^/datum = 7.9 (cf. Table 3). While 7.9 
is clearly unacceptable for any meaningful application, a x^/datum of 1.7 may be 
sufficient for most purposes. 

I turn now to the pp data, Table 4. Typically, x^ for pp data arc larger than for 
np because of the higher precision of pp data. Thus, the Argonne Vis produces a 
X^ /datum =1.4 for the world pp data below 290 MeV and the best Idaho N^LO 
pp potential obtains 1.5. The fit by the best Bochum/ Juelich N^LO pp potential 
results in a x^/datum = 2.9 and the worst is 22.3. 

4. Limitations of xPT 

Since xPT is a low-momentum expansion, we have to expect limitations concern- 
ing its applicability. This is demonstrated in Fig. 2, where phase shift predictions 
by various NN potentials are shown up to 1000 MeV lab. energy for the incident 
nucleon. The figure includes one representative of the family of the high-precision 
NN potentials, namely, the CD-Bonn potential* (solid line), which obviously pre- 
dicts the phase shifts correctly up to the highest energies shown, even though it was 
adjusted only up to 350 MeV. The same is true for other high-precision potentials, 
like Argonne Vig^* and the Nijmegen potentials.^ On the other hand, the chiral NN 
potentials at order N^LO (dashed line Ref. 16 and dotted line Ref. 25) do not make 
any reasonable predictions beyond about 300 MeV lab. energy. This is, of course, 
not unexpected since xPT applies only for momenta Q <C A^ « 1 GeV, which 
is enforced by the regulator Eq. (1) where a typical choice for the cutoff param- 



Table 4. x^/datum for the reproduction of the 1999 pp database by various pp potentials. 
Numbers in parentheses denote cutoff parameters in units of MeV. 







Idaho 


Bochum/ Juelich 


Argonne 


Bin (MeV) 


# of pp 


N^LQi^ 


N^LO^s 






data 


(500-600) 


(600/700-450/500) 




0-100 


795 


1.0-1.7 


1.0-3.8 


1.0 


100-190 


411 


1.5-1.9 


3.5-11.6 


1.3 


190-290 


851 


1.9-2.7 


4.3-44.4 


1.8 


0-290 


2057 


1.5-2.1 


2.9-22.3 


1.4 



February 2, 2008 4:9 



WSPC - Proceedings Trim Size: 9.75in x 6. Sin exoct07 



7 




250 500 750 1000 250 500 750 1000 

Lab. Energy (MeV) Lab. Energy (MeV) 



Fig. 2. np phase shifts of the ^Sq and ^Pq partial waves for lab. energies up to 1000 MeV. The 
solid curve shows the phase shifts predicted by the CD-Bonn potential.* Note that this curve is 
hardly visible because it agrees with the data and, thus, is buried under the symbols representing 
the data. The dashed and the dotted lines are the predictions by the N^LO chiral potentials 
constructed by the Idaiio^® and the Bochum/Juelich^^ groups, respectively. Solid dots represent 
the Nijmegen multienergy np phase shift analysis^® and open circles the GWU/VPI single-energy 
np analysis SMQQ.^'^ 

eter is A « 500 McV. Thus, chiral potentials arc reliable only for CMS momenta 
p,p' < 2.2 fm^^. A Fermi momentum kp w 2.2 fm^^ is equivalent to a nuclear mat- 
ter density p w 4po where po denotes nuclear matter saturation density. Nuclear 
matter calculations in which chiral two-nucleon potentials are applied can be found 
in Rcfs. 28,29. Note, however, that for "complete" calculations also the chiral 3NF 
must be included. In any case, based upon the above arguments, one may trust the 
XPT approach up to densities around 4po- In contrast, relativistic meson theory can 
be trusted to very high momenta (cf. the CD-Bonn curve in Fig. 2) and densities 
equivalent to those high momenta. 

5. Conclusions 

The theory of nuclear forces has made great strides since the turn of the millennium. 
Nucleon-nucleon potentials have been developed that are based on proper theory 
(EFT for low-energy QCD) and are of high-precision, at the same time. Moreover, 

the theory generates two- and many-body forces on an equal footing and provides 
a theoretical explanation for the empirically known fact that 2NF ^ 3NF ^ 4NF. 

At N^LO,^^'^* the accuracy can be achieved that is necessary and sufficient 
for microscopic nuclear structure. First calculations applying the N'^LO NN poten- 
tial^^ in the (no-core) shell model, ■'""''■^ the coupled cluster formalism, and the 
unitary-model-operator approach'^^ have produced promising results. 

The 3NF at NNLO is known^°'^^ and has been applied in few-nucleon reac- 
tions^^''*"'^^ as well as the structure of light nuclei. However, the famous 'Ay 
puzzle' of nucleon-deuteron scattering is not resolved by the 3NF at NNLO. Thus, 
the most important outstanding issue is the 3NF at N^LO, which is under construc- 
tion. 
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